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Abstract 

The concept of p-orthogonality (1 < p < n) between n-particle states is introduced. 
It generalizes common orthogonality, which is equivalent to ra-orthogonality, and 
strong orthogonality between fermionic states, which is equivalent to 1-orthogonality. 
Within the class of non p-orthogonal states a finer measure of non p-orthogonality is 
provided by Araki's angles between p-internal spaces. The p-orthogonality concept 
is a geometric measure of indistinguishability that is independent of the representa- 
tion chosen for the quantum states. It induces a new hierarchy of approximations for 
group function methods. The simplifications that occur in the calculation of matrix 
elements between p-orthogonal group functions are presented. 
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1 Introduction 



Recent studies have aimed at denning a geometric measure of entanglement (see [1] and 
therein). Entanglement is also related to the Von Neumann entropy of reduced density 
operators, however, in the case of quantum systems made of identical particles, it has 
proved important to take apart the uncertainty due to the indistinguishability of the 
identical particles from that due to entanglement [2] . In this work we provide a geometric 
measure of indistinguishability. 

Indistinguishability of identical particles is related to orthogonality properties of Hilbert 
subspaces. It is common knowledge that sets of identical fermions can be considered as 
distinguishable, when their respective wave functions (or density operators) are built 
from one-particle functions belonging to orthogonal Hilbert spaces [3]. In other words, 
when any one-particle state of a set of particles does not overlap with any one-particle 
state of another set, then antisymmetrizing or not antisymmetrizing the tensor product 
of the wave functions of the two sets give the same physical predictions. Such sets of 
fermions are said strongly orthogonal to each other [4,5]. However, as far as we are aware, 
when some one-particle states of the two sets do have non-zero overlap, so that the sets 
become indistinguishable, there is no measure to quantify to which extend the particles 
of both sets are actually mixed. 

More specifically, let, a\, . . . , a 2n be 2n creation operators of orthonormal one-particle (ei- 
ther boson or fermion) states. The n-particle states, a{ ■ ■ ■ aJjO), and, a\ • • ■ a n -i a n+i\tyi 
are orthogonal. So are, a| • - • a^|0), and, a„ +1 • • • a^jO). Intuitively the latter pair is 
"more" orthogonal than the former. In fact, it is "strongly" orthogonal [4,5]. Between 
these two extreme cases, there are intermediate cases, like for example the pairs, a\ • • ■ aJjO), 
and, a\ ■ ■ • a) n _ p a n+1 ■ ■ -a^JO), which are orthogonal but not strongly orthogonal. The 
aim of the present article is to introduce a graded orthogonality concept which discrimi- 
nates between all these cases. Our geometric concept is well-defined (i.e. independent of 
the arbitrarily chosen representation of quantum states) for general, multiconfigurational 
wave functions of possibly different particle numbers as well as for mixed (i.e ensemble) 
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states. 



The article is organised as follow, we first recall the definition of the p-internal space 
of an n-particle quantum state, then we define the concept of p-orthogonality and the 
Araki's angles between the p-internal spaces, finally we show the usefulness of these 
concepts to simplify the calculation of matrix elements appearing in a class of general 
approximation methods for solving the Schrodinger equation of n identical particles. 
Note that throughout the article, the emphasis will be put on fermions, because p- 
orthogonality will be a priori more useful for sets of particles obeying the Pauli principle, 
than for bosonic particles, whose states tend to degeneracy rather than to orthogonality. 

2 p-internal space of an n-particle state 

Let H denotes the one-particle Hilbert space and A n H, (respectively, V n 7i), the Hilbert 
space of antisymmetric, (respectively, symmetric), n-particle states built from H. Let 
^ G A n 7i (respectively, \1/ G V n 7Y) be a normalized n-fermion (respectively, n-boson) 
wave function. Its reduced density operator, D^, acts on a wave function, $ G ATi 
(respectively, $ G V7i) in the following way [7], 

= ^ $ ^ ^ ; (1) 

where (resp. <—>) denotes the right (resp. left) interior product. 

We recall that the interior products for fermions are defined by conjugation with respect 
to the Grassmann product: G A q ~ p H, ^ G A P H, $ G A q H, 

(e|* = a e|$), (2) 

(0|$^ = (9A*]$). (3) 

Similarly, for bosonic states the interior products are conjugated to the symmetrical 
product V. 
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The interior products are equivalent to "annihilation" in the second quantization lan- 
guage; using this formalism £>^($) would be written: Dy\$) = ^(Ol^O^^^O) |0), 



where 0, (bold symbol) denote the annihilation operators associated to 9; ^ , Q\ 
denote their conjugate creation operators. The reduced density operator preserves the 
number of particles, that is to say, A P H, (respectively, \/ p H), is stable under D^. The 
restriction of to the p-particle space, D^,, is the so called "p-order reduced density 
operators" , (the action of D\ can be extended to the whole of AH, (respectively, VH) 
by -D^,($) = if $ G A q H (respectively, VH) and q ^ p, then decomposes as a 
direct sum Dq, = Y, P >o D^,). 

We call "p-internal space" the sum of the eigenspaces of the p-order reduced density 
operator associated to non zero eigenvalues. The p-particle functions of this space 
are called "p-internal functions". An alternative definition of the p-internal space, of 
^ E A n H, denoted l p [^], is: 



that is to say, X P [\I/] is the vector space obtained by annihilating a (n— p)-fermion function 
in ^ in all possible manners. A similar definition holds for bosons. 

Examples from electronic structure theory: The 1-internal space, or simply "internal 
space", X 1 ^], is the space spanned by the, so-called, occupied, natural spinorbitals, in 
quantum chemistry. The 2-internal space, X 2 [^], is the space spanned by the occupied, 
natural geminals. The n-internal space is the one-dimensional vector space spanned by 
the wave function \E'. The p-internal space of a single configuration function (Slater 
determinant) built over a set of n orthogonal spinorbitals, \1/ := 0i A . . . A (j> n , is the 



(^J -dimensional vector space spanned by the p-particle functions, 0^ A . . . A 0« p , built 
over p spinorbitals of ^. 



e 



Zfftf] : = {$ e A P H, 3Q e A n ~ p H, tt ^ V = $}, 



(4) 
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The definition extends to ensemble states described by a general density operators, D n , 
that is, by a convex combination of pure states density operators: 

D n := c i D li with a > and ^ q = I. (5) 

i i 

In such a case, the p-order reduced density operator to consider is, simply, 

= (6) 

i 

It is easy to see that the p-internal space of D n is the sum (not necessarily direct) of its 
pure states p-internal space, 

i p [D n ] = j2i p m- (7) 

i 

The orthogonal complement of the p-internal space, that is the kernel of D p , is called 
the p-external space, £ p [D n ] := X p [D n ] ± , and satisfies, 

£P[D n ] = r\£r[*i]. (8) 

i 

3 p-orthogonality 

3. 1 Definition 

Let G A ni 7i and \& 2 £ A n2 TC be respectively a n\- and a n 2 -fermion wave function. 

We will say that and ^2 are p-orthogonal (for 1 < p < inf(ni,n2)) if and only if 
their p-internal spaces are orthogonal, 

J p [^] JL 1 P [V 2 }. (9) 
A similar definition holds for bosonic states, and extends to ensemble states, either 
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bosonic or fermionic, by considering the orthogonality of the p-internal space of their 
associated density operators, D™ 1 and D 2 2 . 

We see immediately that if n\ = n 2 = n, n- orthogonality is the usual orthogonality 
between wave functions. In the case of ensemble states, it means that any wave function 
associated to a pure state in the convex combination of one density operator is orthogonal 
to any wave function associated to a pure state in the convex combination of the other 
density operator. 

At the other end, 1-orthogonality between and ^2 amounts to strong orthogonality, 
usually defined by, 

J dTi*i(Ti,T 2 ,...,T n J* 2 (Ti,T2,...,< 2 ) = Vt 2 , . . . , T ni , T 2 , . . . , T^. (10) 

This can be rewritten, using Dirac distributions centered on the Fermion variables, as 
the nullity of the kernel, 

(5 T2 A ... A 5 Tni <- A ... A 5 <2 <- # 2 > = 0. (11) 



or, by changing to a basis set representation {<pi}i in the rigged Hilbert space [8], as, 

(<f>i 2 A ... A (f) ini ^i\4>i' 2 A ... A i / >2 ^* 2 ) = Vi 2 ,...,i ni ,i 2 ,...,C (12) 



Since, the (ni — l)-particle functions, (f>i 2 A . . . A <pi n , span all of, A^-^Ti, and the 
(n 2 — l)-particle functions, 0^ A ... A 0^ , span all of, A^ n2 ~^H : the latter equation is 
equivalent to orthogonality between any pair of 1-internal functions, that is to say to 
1-orthogonality. 

Remark: The present definition of "strong orthogonality" as the orthogonality of the 
one-internal spaces, and, another characterization of the one-internal space of a function 
^ G A n H as the smallest Hilbert space, T such that \1/ G l\ n T [9] , make obvious the sep- 
arability property of strongly orthogonal electron pairs [6], or more generally, of strongly 
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orthogonal electron groups (see also the definition of "non-overlapping subsystems in 
[10]). 



3.2 Graded orthogonality 

An important property to notice is that p-orthogonality implies g-orthogonality for all 
q > p. Its proof relies essentially upon the following lemna: 

Lemna 1: If $ G T p+1 [D n ] then V0 G^, (^$)6 l p [D n ]. 

Proof: Consider first the case of fermionic pure states. Let \l/ G A n 7i be an n-particle 
wave function, $ G J p+1 [^] and G TL. By Eq. (4), there exists Q G A n ~ v ~ x, H such that 
il ^ # = $. So, $ = (ft *) = (n A 0) V, where (ft A 0) G A n ~ p 7Y. 
This means that (0 $) G J 7 ^] according to Eq. (4), and proves the proposition 
for fermionic pure states. The demonstration is the same for bosonic pure states with V 
instead of A. 

Now consider a mixed state operator D n as in Eq. (5), $ G T p+1 [D n ] and G 7i. By Eq. 
(7), there exist $,'s such that $ = £ $j and G J p+1 [^j] for all i By (anti) linearity 

i 

of the interior product, $ = (j2 ^^j — Z) (0 ^ But we have just shown 
that (0 <— > $j) G X p [^j] for all i, which proves the property for mixed states according 
to Eq. (7). 

p-orthogonality is a graded property in the sense that: 

Proposition 1: If two states represented by the density operators D™ 1 and D% 2 , (or by 
the wave functions \&i and * 2 for pure states, with D™' = |\l>j)(\l/j|), are p-orthogonal 
then they are a fortiori g-orthogonal for all q such that, m/(ni,n 2 ) > q > p. 
Proof: 

Let 1 < p < rii < n 2 be three integers. Let {0i}« be an orthonormal basis set of 
H. Consider the fermionic case, and note that the "particle number" operator, N : = 
J2i4>i A (0j ^ •), acts 011 the g-fermion space, A q H, as, q.ld^u-, {Id/\iH denotes the 
identity on A q H). For all T 1 G 1 P+1 [D^}, T 2 G X P+1 [D^ 2 ], and & G W, we have 
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((Pi <- ri|0i r 2 ) = since by lemna 1, <- r\) G 3?[D?], (0, T 2 ) G T?[D?\, 
and by hypothesis -D™ 1 and _D 2 2 are p-orthogonal. Therefore, = J2i(4>i ^ ^i\4>i ^ 

r 2 ) = E,(0, a (^roira) = (ivr 1 |r 2 > = ( P + lX^r,). So, vr\ g t*\d?\, 

Vr 2 G IP+^DJ 2 ], (ri|r 2 ) = 0, which proves that, D™ 1 and D 2 2 are (p + l)-orthogonal. 
The proof works for bosons, if A is replaced by V. By induction, the result holds for all 
q such that, inf(ni,n 2 ) > q > p. 

So, p-orthogonality provides us with a graded orthogonality concept for states of identical 
particles, and the traditional term of "strong orthogonality" attached to 1-orthogonality 
is justified in the sense that it implies p-orthogonality for all p. 

Example 1: For integers, n > p > 0, the pairs, \Pi := 0i A . . . A n and \& 2 := <pi A . . . A 
4> n -p A <pn+i • • • A <pn+p, (equivalent to those denoted with second quantization operators 
in the introduction), are (n — p + l)-orthogonal but not (n — p)-orthogonal since for 
$1 := 4>n- P +i A ... A n and $ 2 := n+1 A ... A (j) n+p , ($i #i|$ 2 W 2 ) = (0i A . . . A 
<f>n- P \<f>i A ... A 4> n -p) = 1 is non zero, although ($; vjr) G J n_p [\I>i], for i G {1, 2}, by 
definition. 

Example 2: Let (0j)j=i,...,8 be 8 orthogonal spinorbitals. The functions := 0i A 2 A 
03 + 04 A 05 A 6 and ^ 2 := 0! A0 7 + 2 A0 8 are 2-orthogonal (it is impossible to obtain \1' 2 
by annihilating a spinorbital in vE^) but not 1-orthogonal since both 0! and 2 belongs 
to their one-internal space. 

3.3 Araki angles 

Within a given "graduation", e.g. the set of functions which are (p+ l)-orthogonal but 
not p-orthogonal for some p, a finer measure of non p-orthogonality is given by the Araki 
angles between the p-internal spaces. The Araki angles between the spin a- and the spin 
/3-part of the one-internal spaces have already been introduced by the present author to 
study spin contamination in spin-unrestricted wave functions [11]. The cosines of these 
angles are the overlaps between biorthogonal functions. 
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Consider the p- internal spaces P^P™ 1 ] and Z V \D\ 2 \ of two density operators (or wave 
functions in case of pure states), with n\ < n 2 - Let us set E := JPfP™ 1 ] + J p [£>2 2 ] an d 
denote Pj (j £ {1, 2}) the orthogonal projector on J P [P™ 3 ] in E. The construction is the 
same as that of [11]. We define the operators, u COSQ pv and "SINQ?", 

COSQ p := | Pi + P 2 - Id E \, SINQP := |P X - P 2 |, (13) 

which satisfy, 

{COSQ p ) 2 + {SINQ P ) 2 = Id E . (14) 

{COSQp) 2 is a Hermitian, positive operator, whose eigenvalues are in the interval [0, 1]. 
One can associate to each eigenvalue, \ P , an angle by, 

0f = arccos( y f>? i ). (15) 

The eigenspaces of {COSQ p ) 2 , that we write V e v (rather than V x p) decomposes E into a 
direct sum of orthogonal vector subspaces, 

£:=0V>. (16) 

The Araki angle operator, Q p , is defined on E as 

6*:=5: flf.JV (17) 

i i 

where Pv gP is the orthogonal projector on V 8 p. The remarkable property of the decompo- 
sition (16) is that it "respects" the structure of the p-internal spaces J^P™ 1 ] and J^fP^ 2 ], 
in the sense that, for j e {1,2}, 

lP[D^] = ^lP[D^]f]V eP . (18) 

e p 

Setting ZP\D]%v := lP[D r - 3 ]f]V e P we obviously have that 1 p [D^} 9 p is orthogonal to 
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T p [D% 2 ] e P if i ^ j, and if i = j, any pair of function $1 G J^D™ 1 ]^ , $ 2 G X p [£>2 2 ]^ can 
be thought geometrically as making an angle Of . 

Particular cases: 

If there is only one eigenvalue, \\ = 0, hence Q\ = \, the states are in fact p-orthogonal. 
If the eigenvalue, Af = 1, hence 9f = 0, is present and its multiplicity equal to ni, then 
the p-internal space X p [Di l ] is a vector subspace of J^Z^ 2 ]- 

Between these extreme cases, the Araki angles provide us with a quantitative mean to 
assess departure from p-orthogonality. 



4 Application to group functions 



When the respective states of two groups of identical particles are p-orthogonal, at most 
p—1 particles are possibly "overlapping" over the two groups, in the sense that the overlap 
of any p-particle state occupied in the ni-particle state of the first group with any p- 
particle state occupied in the n2-particle state of the other group is zero. In particular, 
for p = 1, no particle overlaps and the two groups are distinguishable. 

Let us emphasize that this notion of distinguishability does not necessarily imply the 
localization of the two groups of particles in two non-intersecting regions of real space. 
It has only to do with the orthogonality of abstract Hilbert spaces. 

A direct consequence of p-orthogonality is the cancellation of matrix elements between 
Hermitian operators that only couple a limited number of particles: 
Proposition 2: Two p-orthogonal, n-particle states cannot be coupled through a q- 
particle interaction operator, V q , if q < n — p. 

Proof: A g-particle operator, V q , is an operator that can be expressed in the second 
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quantization formalism as, 

V*= £ \i,jai...a\a jq ...a h , (19) 

J- = (h,-,jq) 

with A/j = A} 7 . Let and \l/2 be two p-orthogonal, n-particle wave functions. By 
linearity, 

(tfxll^tta) = E A/,j (^ila^ . ..a} a,- . . .a n \V 2 ) = E A/,j (a*, . ..a^a^ . ..a h V 2 ). By 

i, j q i, j 

definition, V(ii, . . . , i,), (ji, • ■ • , j q ), Oi, • • ■ e a jq . . . e J"-«[* 2 ]. 

By hypothesis, p < (n — q), so Proposition 1 shows that all these pairs of (n — g)-particle 
wave functions are orthogonal: 

V(ii,...,g,0'!,...,i 9 ), (a iq ,..., 0il ^i| ^-,,...,0^*2) = 0, hence (*i|y 9 |* 2 ) = 0. 

In quantum chemistry, general antisymmetric product function methods [4,5,12,13,14,15,16,17,18,19,20,21,1 
optimize n-electron wave functions of the form: 

$ = $1 A ■ ■ ■ A $ r , (20) 

where is an n^-electron function and Ei n i — n - So far, for practical purposes, all these 
approaches (except those of [27,36]), have imposed the constraint that the ^j's have to 
be 1-orthogonal to one another. In the Electronic Mean Field Configuration Interaction 
(EMFCI) approach [36], no orthogonality constraint is a priori imposed on the In 
particular, in the simple case where, for all i, rii = 2, both APSG (Antisymmetrized 
Product of Strongly orthogonal Geminals) [31] and AGP (Antisymmetrized Geminal 
Product) of extreme type [38], $1 A ■ • • A ^1, are considered by the EMFCI variational 
process. 

Therefore, it would be interesting to analyse the optimized EMFCI functions obtained 
for different systems and geometries in terms of their p-orthogonality properties, and 
see for instance, if they are closer to the APSG case (1-orthogonality) or to the AGP of 
extreme type case (non 2-orthogonality with the Araki angle equal to zero for all pairs 
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of two-electron group functions). However, in the present study we will limit ourselves 
to emphasize how enforcing a p-orthogonality constraint between the \IVs, simplifies the 
computation of the Hamiltonian and overlap matrix elements between general antisym- 
metric product functions (Eq.(20)). 

Let us consider another such function, \P' = ^[ A • • • A ^' r , (with n\ = rij), and define the 
notation, 

*;:=*iA-"A* j _iA* i+ iA---A* ri and so on, 

to denote that one or more specified factors have been taken out of a product function. 
So, for example = ty[ A ^j. It can be shown, using the Hopf algebra tools of [36] , 
that, 

(*i A tyjl^i A • • • A * r ) = 

E P\n\, ni -\n..,\nn r -\ir\ a • • • a (* r )p.) (^1(^1)71 a • • • a (* r ) Jr ) 

J 1 , ...,/-• 

Vj |7^|e{0,...,n,}, 

(21) 

where 1 , 

EEl^l-K-l^l) 

P|/l|,n 1 -|/l|,...,|/1,n r -|/''| = (-l)''= afc=1 , (23) 



and where, for any p-particle wave function, $ := E \k V'fci A • • • A (X 

A":=(fcl<...<fcp) 

runs over ordered sequences of positive integers, (^)j denotes a one-particle basis set), 



1 this formula is a particular case of a formula given in [36,40] with an error on the summation 
bounds. A correct version is, 

E E « 

Pn\,... K ,...,nl,...,nl = (-1)^' **>^ (22) 
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and any ordered sequence of length m G {0, . . .p}, I :— (1 < i\ < . . . < i m < p), the 
following compact notation is extensively used, 

• • • ($) 7 • • • ($)/•• • := Pl j J2 X k ■■■ A • • • A Vfe im ) 

K:=(k 1 <--<k p ) 

■ ■ ■ A • • • A ^ ) • • • (24) 

with 7 := (1 < ii < . . . < i p - m < p), complement of / in {1 < 2 < • • • < p}, pjj is the 
sign of the permutation reordering the concatenated sequence ////in increasing order; 
if the length, |/|, of J is then, by convention, ($)/ := ($)0 = 1, and P0,(i<...< p ) = 1; 
note that, ($)(i<... <p ) = $. 

If we assume that the group-1 (called the group of active electrons in the EMFCI method) 
wave function, ^f[, is g-orthogonal to the product of the wave functions of the other 
groups (called spectator groups in the EMFCI method), = $2 A • • • A \l/ r , Eq.(21) 
becomes, 

(*i A ^jl^i A • • • Atf r ) = 

E p\nn,-\n..,\nn r -\ir\ WK^iV a • • • a (* r )r-X*'il(*i)/i a • • • a (* r 

i 1 ,...,r 

^^{m-q+l,...^!}, 

Vj>i \P\e{o,..., nj }, 

£ |/ J '|=ni 
i=i 

(25) 

that is to say, the summation on the ordered sequences, I 1 , is limited to those with length 
strictly more than n\ — q. Without this restriction, the number of /^sequences would 
^ e (rTi) = S (T) (n7-i) ' wnereas with the g-orthogonality restriction, it falls down to 
£ f™ 1 ) (n~-i) • In tlie limit case of 1-orthogonality, only the sequence I 1 = (1 < 2 < 

i=m— q+l 1 ' *' 

• • ■ < rix) remains, and Eq.(25) simplifies to, 

(*i A^'jl^i A---A* r ) = (* 1 |*i)(* / i |* 2 A---A* r ). (26) 
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At the other end, enforcing ni-orthogonality between active and spectator groups, that 
is the weakest g-orthogonality constraint, rules out only the case J 1 = 0. However, this 



of active electrons, ni, is small with respect to the total number of electrons, n. 

Consider now the matrix elements between general antisymmetric product functions, of 
an operator, H, whose action on n-electron wave functions is induced by a s-particle 
operator, h, (with s < n). Typically, h will be a Coulombian Hamiltonian, so that s — 2. 
Its induced action on the n-electron wave function of Eq.(20) can be expressed using 
Hopf algebra techniques as, 

if [tf i A • • • A * r ] = 

E P\JW,...,\Jr\,\Jr\(*l)jl A • • • A (V r )jr A h J! A • • • A (* r ) Jr] . (27) 

j\...,J r 

\/j \j*\e{o,..., nj }, 

So, a matrix element, A ^j|if[^i A ■■• A ^ r ]), is a sum of terms of the form, 
(*i A ^K^Oji A • • • A (^ r )jr A h [(*i)ji A • ■ • A (* r )jr]), similar to Eq.(21) but with, 
in the ket, (r + 1) groups of | J 1 !, . . . , | J r \, s particles respectively, instead of r groups of 
ni, . . . ,n r particles. In particular, | J x | can be less than ni, the number of particles in 
However, essentially the same development can be carried out, the g-orthogonality 
constraint limiting the summation for each term, 
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(*i A ^K^Oji A • • • A (* r )jr A h [(^Oji A • • • A (* r )>]) = 

12 P\n\,\n\,...,\ir+i\,\ir+i\ 

|/1|e{ ;;- / ;; i .. i|J1|} • m WiM* a ... a ((*,)^ a k*^) a ... a (^),>]) jr+1 > 

l<j<r, |P'|6{0,...,|^'|} 

|r-+i|6{o,... )a }, • ((*i)ji)ji A ••• A ((%)jr) Ir A (/i [(^i)ji) A 

r + 1 



A(^)jr]j /r+1 ). 

(28) 



5 Conclusion and prospects 



We have defined the geometric concept of p-orthogonality between quantum states of 
sets of identical particles. This concept provides us with a graded measure of indistin- 
guishability in the sense that two sets of identical particles that are g-orthogonal can be 
seen as "more indistinguishable" than two sets that are p-orthogonal if q > p, because a 
larger subset of particles can possibly share i.e. occupy a substate, common i.e. internal 
to the quantum states of both sets. 

A classical anology can be attempted with the case of two groups of billiard balls of the 
same color. When p — 1, no particle is mixed and the two sets of particles are in fact 
distinguishable like two sets of balls localized in distinct areas of a billiard table. Pushing 
this classical picture one step beyond for p > 1, the two sets of balls would be connected 
but at most p — 1 balls of one set would be in contact with at most p—1 balls of the other 
set. So, the smaller p, the narrower the bridge between the two sets of balls would be. 
Assuming that the group of origin of the balls making up the bridge is unknown, these 
balls would be the analogues of the genuinely indistinguishable particles which belong 
partially to both sets. 

Let us emphasize that this classical picture should not be carried too far, for, in particu- 
lar, our notion of distinguishability does not imply the localization of the two groups of 
particles in two non-overlapping regions of real space. It has only to do with the orthogo- 
nality of Hilbert spaces called the p-internal spaces of the quantum states, p-orthogonality 
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can be seen as a mathematical and quantum mechanical rigorous formalisation of this 
classical image. 

p-orthogonality can be used to remove some arbitrariness in the choice of a representation 
for a quantum system in the same manner as localization criteria do. For example, con- 
sider n pairs of spin-^ fermions whose state can be represented by a Slater determinant, 
^ := 0i A 0i A • • ■ A (fi n A (j) n , where <f>i, . . . , (fi n , are orthonormal one-fermion functions of 
spin ,2-component |, and (fii, . . . , (fi n , their counterparts with spin z-component equal to 
— | . Such a wave function is invariant within a phase factor under an unitary transfor- 
mation, u, of the one-particle functions, <f>i, . . . , <f> n . There are various techniques [42,43] 
that exploit this freedom to reexpress a wave function with a new set of one-particle 
functions, ifij = u(<j>j), localized in real space, and such that \I/ has still the form of a 
Slater determinant, \I/ := ifti A'lfii A • • • Aip n Aip n . However, this only provides a constraint 
on one-particle states and there is still more freedom available. For example, \I/ can be 
re-expressed as an AGP of extreme type with the same localized one-particle functions, 

= g A g A ■ ■ ■ A g , (29) 

nf actors 

where g = (n\)~n {jpi A ty\ + • — h ip n A ■?/>„). If we set, gi — tfii A ifii for all i, we also have 

* = 9i A <? 2 A • • • A g n . (30) 

Imposing 1-orthogonality or even 2-orthogonality between the two-fermion functions 
appearing in Eqs.(29) and (30) can discriminate between these two equivalent writings. 

The graded structure of p-orthogonality constraints naturally leads one to consider a 
corresponding hierarchy of approximations for methods based on general antisymmetric 
product functions. In this work, we have exhibited the link between p-orthogonality and 
the combinatorics involved in the calculation of the matrix elements of particle-number- 
preserving observables. In the frame of the EMFCI method, we have shown that even 
the weakest p-orthogonality constraint between an active group of particles and the rest 
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of spectator particles can be effective in limiting the computational effort required for 
the calculation of matrix elements. We will report shortly on the accuracy of EMFCI 
wave functions constrained by p-orthogonality, for increasing value of p on a benchmark 
of molecular systems. 
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